Universal fluctuations and extreme value statistics 
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We study the effect of long range algebraic correlations on extreme value statistics and demonstrate 
that correlations can produce a limit distribution which is indistinguishable from the ubiquitous 
Bramwell-Holdsworth-Pinton distribution. We also consider the square-width fluctuations of the 
avalanche signal. We find, as recently predicted by T. Antal, M. Droz G. Gyorgyi and Z. Racz 
for logarithmic correlated 1/f signals, that these fluctuations follow the Fisher-Tippett-Gumbel 
distribution from uncorrelated extreme value statistics. 
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I. INTRODUCTION 

Three years ago Bramwell, Holdsworth and Pinton 
(BHP) published the remarkable discovery that the same 
functional form that describes the fluctuation spectrum 
of the energy injected in an experiment on turbulence 
also describes the fluctuations in the magnetization of 
the finite size two dimensional XY equilibrium model in 
the critical region below the Kosterlitz-Thouless transi- 
tion temperature Q . Since three dimensional turbulence 
and two dimensional magnetic equilibrium system appear 
to have very little in common, BHP made the reason- 
able suggestion that the origin of the identical functional 
form for the fluctuation spectra should be sought in the 
one thing the two systems appear to share, namely, scale 
invariance. This suggestion was supported by the subse- 
quent finding that a long list of scale invariant (or nearly 
scale invariant) non-equilibrium as well as equilibrium 
systems exhibit the same BHP form for the fluctuation 
spectrum for certain quantities |^| . 

Nevertheless, not all critical systems fluctuates accord- 
ing to the BHP spectrum. This was made very explicit 
by Aji and Goldenfeld ||, who proffered the interest- 
ing suggestion that the reason the two dimensional (2d) 
XY-model and the driven turbulence experiment exhibit 
the same fluctuation spectrum could be that the 2d XY- 
model is the effective model for the turbulence experi- 
ment. Even if this is correct we still lack an explanation 
of why so many disparate systems Q do exhibit BHP 
fluctuations. 

The BHP functional form is similar to one of the 
asymptotic forms for extreme value statistics: the asymp- 
tote first discussed by Fisher and Tippett Q and often 
referred to as the Gumbel distribution || . The BHP dis- 
tribution is, however, not identical to the Fisher-Tippett- 
Gumbel (FTG) asymptote. Consider T independent 
and identically distributed stochastic variables. Under 
certain conditions (essentially exponential tail) the fc-th 
largest of the T variables will be distributed according 
to the FTG asymptote with k entering as a parameter. 
The BHP form can be thought of as corresponding to the 
somewhat uninterpretable case of k = n/2 

As already alluded to in Ref. Q the deviation between 
the BHP and the FTG form may be related to correla- 



tions. Our main aim in the present paper is to study this 
point in detail. We do that by simulating the so called 
Sneppen depinning model [Q . The power spectrum of the 
Sneppen model behaves like 1/ '/^ with (3 ~ 0.5 for low 
frequencies corresponding to a very slow algebraic decay 
of the autocorrelation function (see e.g. ^]). We use the 
Sneppen model in the next section to demonstrate that 
extreme value statistics of T strongly correlated expo- 
nentially distributed variables may follow the BHP form. 
Remarkably, Antal et al. [|| demonstrated analytically 
that, at least for a certain class of 1/f signals (peri- 
odic signals), the width-square fluctuations (W2) of the 
signal follow the FTG distribution from extreme value 
statistics, though it is not clear why this should be the 
case. Inspired by this finding we study in Section III the 
width fluctuations of the avalanche signal in the Snep- 
pen model. We find, contrary to the 1/f case studied 
by Antal et al., that the probability density function for 
W2 in the Sneppen model is very well represented by the 
FTG distribution even for non-periodic signals. Section 
IV contains a discussion and our conclusions. 



II. EXTREME VALUE STATISTICS AND THE 
SNEPPEN MODEL 

To investigate the relationship between extreme value 
statistics of correlated variables and the BHP probability 
density we consider now the simple 1+1 dimensional de- 
pinning model introduced by Sneppen [Q. The model is 
imagined to represent a one dimensional elastic interface 
moving transverse while acted upon by a set of random 
pinning forces, See Fig. la. 

The model, which is discrete and very schematic, con- 
sists of L sites in the x-direction and infinitely many sites 
in the y-direction. Each square on this semi-infinite lat- 
tice is assigned a random number, the pinning force, uni- 
formly distributed on the interval [0, 1[. The interface is 
represented by the set {(x, y)\y = h(x,t)}, where h(x,t) 
denotes the "hight" of the interface above the x-axis at 
time t at site x along the x-axis. The initial configura- 
tion is h(x, 0) = for V x. In each time step the interface 
site, x s with the smallest pinning force is located and the 
interface at this location is moved one step ahead, i.e., 
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h(x s ,t) t— > h(x s . t) + l. This may cause the neighbouring 
slope \h(x s ,t) - h(x s - l,t)| (or \h(x s + l,t) - h(x s ,t)\) 
to exceed 1, in which case the interface at site x s — 1 is 
moved one step ahead, i.e. h(x s — 1, t) t— > h(x s — l,t) + l 
(and similar for the site x s + 1 if needed). The update 
of the nearest neighbour sites of x s may cause the slope 
on the next nearest sites to exceed 1. In which case the 
interface is moved ahead on these sites. This procedure 
is repeated until all slopes satisfy h(x, t) — h(x—l, t)\ < 1 
once again. The sequence of operations needed to make 
all slopes smaller than or equal to 1 after the update of 
site x s is denoted one time step. The number of sites up- 
dated during one time step is called a micro-avalanchc. 
The number of sites being updated during the time step 
t is called the size of the avalanche s(t) and per defini- 
tion the duration of each of these avalanches is the same, 
namely one time step. 
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FIG. 1. a) Snapshot of the sneppen interface h(x) at 
two times 10000 timesteps apart, system size L — 500. b) 
The PDF of the micro- avalanches in the sneppen model, 
system size L = 500. c) The power spectrum of mi- 
cro-avalanches for 4 different system sizes. From bottom to 
top L = 50, 100, 500, 1000. 

The probability density function (PDF) for the size of 
the micro-avalanches is close to an exponential as seen in 
Fig. lb and first reported in Q. The individual micro- 
avalanches are strongly correlated in time. This is clearly 



seen from the power spectrum S(f) of the temporal sig- 
nal s(t), see Fig. lc. The low frequency behaviour of 
the power spectrum is approximately S(f) ~ 1/V7> in- 
dicative of slow algebraic decay of the auto-correlation 
function of the signal s(t) ~ l/V* ( se e e.g. M). 

Let us now imagine that the activity of the model is 
monitored by some devise which has only a limited res- 
olution. This might be modeled by assuming that the 
measured signal St(t), rather than consisting of the mi- 
croscopic instantaneous activity s(t), is given by the sum 
of activities within a time window of a certain size T: 

s T (.T) = jr s (t). (i) 

t—T 

This signal was among the set of quantities in Ref. Q 
found to follow the BHP distribution for a certain range 
of window sizes T. That is, if T is too small the PDF for 
St{t) is close to an exponential distribution. The central 
limit theorem will, however, cause the PDF for St(t) to 
approach a Gaussian when T becomes large enough that 
correlations amongst the s(t) entering the sum in Eq. |l| 
can be neglected. This deviation can be observed in Fig. 
2a. Presumably the PDF for S(T) would remain close to 
the BHP form even for T — > oo if the correlation time of 
the signal was infinite as it would be expected to be in 
the limit of infinite system size, see Fig. lc. 

Since, as seen in Fig. 2b, the functional form of the 
BHP distribution is close to the first FTG asymptote for 
extreme value statistics we proceed to investigate the ex- 
treme statistics of the correlated variables s(t) generated 
by the Sneppen model. For the time windows considered 
for St(t) we define 

M t (t) = max{s(r), s(r + 1), s(r + T)}. (2) 

In Fig. 2b we show the scaled PDF for Mt(t). First 
we compare the distribution of the sum and the max in 
Fig. 2a. We see that the distribution for Mx(t) remain 
very close to the BHP distribution for all considered sam- 
pling sizes T, while the distribution for St(t) gradually 
deviates from the BHP as T is increased. In Fig. 2b we 
demonstrate two interesting points. First that for sys- 
tem size L — 1000 we are unable to make T so large 
that the distribution for Mt(t) deviates form the BHP 
form. No difference between the results for T = 5000 and 
T = 10000 can be detected. The size dependence is, how- 
ever, detectable, as illustrated in the insert to Fig. 2b. 
We show here the PDF for smaller system sizes L = 500 
and L — 50. The spatial extend of the system L = 50 
is now sufficiently small to destroy the long range tem- 
poral correlations as is seen from the fact that the PDF 
for L = 50, T = 500 follows the FTG form for uncorre- 
cted extreme value statistics. We elaborate on the role 
of the correlations in the main frame of Fig. 2b. This 
figure demonstrates that the deviation between the PDF 
for Mt(t) and the FTG form is indeed caused by the 
correlations of the Sneppen model. This conclusion is 
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reached in the following way. We generate T uncorre- 
cted stochastic variables Xi au drawn from the PDF for 
the individual micro-avalanches s(t) (see Fig. lb). The 
difference between 



M uc = max{xi,X2, -,Xt} 



(3) 



and Mt(t) is solely the correlations amongst the primary 
variables from which Mj(t) is generated. Hence, given 
the exponential form of the PDF for the distribution of 
the individual micro-avalanches (see Fig. lb) we expect 
the PDF for the uncorrelated extreme M uc to follow the 
FTG asymptote for large values of T. This is exactly 
what is found in Fig. 2b. 
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FIG. 2. a) The PDF of St(t) for T = 50(+), T = 500(:r) 
and T = 500(*) and the PDF of M T (r) for T = 50 (o), 
T = 500 (□) and T = 5000 (o), the solid line is the 
BHP PDF. b) The PDF of Mt(t) for T = 5000 (o) and 
T = 10000 (A) and the PDF of the decorrelated maximum 
M u c for t = 5000 (*) and T = 10000 (+), the solid line is the 
BHP form and the dashed line is the FTG form. The inset: 
The PDF of M t (t) for system sizes, L = 50, T = 50 (o) and 
T= 500 (•). 



III. THE SQUARE WIDTH AND EXTREME 
VALUE STATISTICS 

In the present section we analyse the statistics of the 
fluctuations in the width of the avalanche signal s(t) gen- 
erated by the Sneppen model and of the 1/ / signal F(t) 
in Eq. o below. Our inspiration is from the beautiful 



paper by Antal, Droz, Gyorgysi and Racz (ADGR) [||. 
These authors demonstrated analytically that at least for 
a certain class oil/ f signals h(t), the width-square 



w 2 (h) = [h(t)-hf 



(4) 



is distributed according to the FTG distribution. In Eq. 
the over-bar denotes the following time average 



(•5) 



t=i 



The result by ADGR is striking since it is unclear how 
or if extreme value statistics is involved in some effective 
way in determening the distribution of W2- 

As a contribution to the investigation of the generality 
of the discovery by ADGR we studied wi for the signal 
s(t). This signal is not 1// but rather and differs 

in this way from the class studied in Ref. Q. As see in 
Fig. 3a W2 for s(t) does follow the FTG functional form 
except for very small sample sizes. This indicates that 
the ADGR result is more general than their calculation 
allows one to conclude. ADGR predicted analytically 
that this should be the case for periodic 1// signals and 
they report that deviations are observed in simulations 
of non-periodic 1// signals. 

To address the question about deviations from the 
FTG form we show in Fig. 3b the fluctuations in w 2 
for a 1// signal of the type considered by Antal et al. 
We generated the signal as the "half integral" of white 
noise M, i.e., the signal is obtained as a convolution be- 
tween white noise and a propagator G(t) which decays 
as an inverse square root: 



F(r) 



lim 

T„->oo 



t=T-T, 



G(r-t) X (t), 



(6) 



where G(x) = l/y/x. We assume to be uncorrelated 
and uniformly distributed on the interval [—1,1] and in 
our simulations we truncate the sum in Eq. (n) at T* = 
2 20 . 

One can calculate the PDF of the Fourier coefficients of 
the signal of length T in Eq. (^|) by performing a discrete 
Fourier analysis (note this amounts to assuming a signal 
of periodicity T) 



FH = 5>(r)e- 



(7) 



T = 



Averaging over the white noise x(i) gives the following 
PDF for the Fourier amplitudes 

rt\fi{ U )\ a =tf) = ^ev(-^I?y (8) 

which demonstrates explicitly that the signal in Eq. (|^) 
is of the type considered by ADGR B . 
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We notice that the simulated distribution for W2 in Fig. 
3b for all sample sizes T deviates from the FTG form pre- 
dicted by Antal et al. This was indeed noticed by ADGR 
(their Fig. 1) and ascribed to the difference between the 
assumed periodic boundary condition needed to perform 
analytical calculation of the PDF of W2 and the non- 
periodic signal simulated. This problem appears not to 
arise in case of the Sneppen model, where Fig. 3a shows 
that the PDF for W2 follows the FTG form as soon as 
T is of order 500 or larger. We recall that the difference 
between the signal in the Sneppen model and the signal 
in Eq. (^) is that the Sneppen signal is algebraically cor- 
related with an autocorrelation function that decays like 
one over the square root of time whereas the 1// signal 
in Eq. (|^) decays even more slowly, namely logarithmi- 
cally. It is likely that this difference in the range of the 
correlations is the cause of the difference between the two 
signals sensitivity to boundary condition. 
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FIG. 3. a) The PDF of the width-square for the sneppen 
model, sample sizes T = 50 (dotted line) and T = 500 (dashed 
line), the solid line is the FTG form, b) shows the PDF of the 
width-square for the 1/f-noise, sample sizes T = 64 (dotted 
line), T = 4096 (dashed line) and T = 32768 (dashed-dotted 
line), the solid line is the FTG form. 



IV. DISCUSSION AND CONCLUSION 

We have studied the effects of correlations on extreme 
value statistics. Our aim has been to relate the Bramwell, 
Heldsworth and Pinton distribution GL^] to correlated 
extreme statistics. We have achieved this successfully for 
one specific system: the Sneppen depinning model. 

We found that the extreme value statistics of the alge- 
braically correlated avalanche signal is described by the 
BHP distribution. Furthermore we studied the width- 
square of the same signal and found this to be dis- 
tributed according to the Fisher-Tippett-Gumbel distri- 
bution for uncorrelated extreme value statistics, though 
no extremes were explicitly involved for the width-square 
signal. This finding generalises a recent result by An- 
tal, Droz, Gydrgyi and Racz ||. The BHP distribu- 
tion and the FTG distribution are of similar functional 
form though they differ significantly for large fluctuations 
away from the mean, see jl2| . 

Extreme value statistics for I// correlated signals were 
studied using renormalisati on g roup (RG) techniques by 
Carpentier and Le Doussal Efj]. They find that the cor- 
relations make the functional form change from the ex- 
ponential decay of the uncorrelated FTG to the form 
yexp(— y) for fluctuations y — x — (x) ^ above the 
average (in the case of Maximum statistics). Note the 
yexp(— y) decay is identical to the asymptotic behaviour 
of the BHP distribution for large deviations above the 
average [pyi^]. Hence, the RG calculation by Carpentier 
and Le Doussal is an analytic indication of a relationship 
between the BHP distribution and extreme value statis- 
tics of logarithmically correlated variables. Our numeri- 
cal study indicates a relationship in the case of algebraic 
correlations, at least in the specific example of the Snep- 
pen model. 

Signals with hierarchical correlations were considered 
in a recent preprint by Dean and Majumdar They 
find that correlations typically, though not always, al- 
ters the super-exponential exp(— exp(y)) found for y < 
in the uncorrelated FTG well as in the BHP 

distribution. This result makes it clear that not all cor- 
relations make extreme value statistics follow the BHP 
distribution. 

We conclude that correlated extreme value statistics 
may produce a limit distribution indistinguishable from 
the Bramwell-Holdsworth-Pinton distribution. More re- 
search is needed in order to determine the generality of 
this result. 
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